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ABSTRACT 


The  pressure  amplitude  and  phase  distribution  along 
the  interface  between  a tapered  fluid  layer  and  an  under- 
lying fast  fluid  bottom  were  investigated  both  theoretically 
and  experimentally. 

Two  different  theoretical  models  were  compared  experi- 
mentally: a simple  model  based  on  a combination  of  normal 

modes  and  ray  theory  and  an  exact  solution  based  on  the 
method  of  images.  The  experiment  was  conducted  at  100  KHz 
with  a wedge  of  silicon  oil  separated  from  a large  tank  of 
fresh  water  by  a thin  mylar  diaphragm.  The  simple  model 
failed  to  predict  adequately  the  pressure  amplitude  and 
phase  along  the  interface.  The  method  of  images  gave 
accurate  predictions. 
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I.  INTRODUCTION 


As  sound  from  the  open  ocean  propagates  up  onto  the 
continental  shelf,  lower  frequencies  become  lost  to  the 
water  column;  energy  is  transmitted  into  the  bottom.  If 
the  bottom  has  a speed  of  sound  in  excess  of  that  of  the 
water,  then  it  is  predicted  [Ref.  1,2]  that  a highly 
collimated  beam  will  form  in  the  bottom. 

For  a channel  of  constant  depth  overlying  fast  bottom, 
normal-mode  theory  shows  that  if  the  frequency  of  the  sound 
is  less  than  the  cutoff  frequency  of  the  mode  energy  can 
be  transferred  into  the  bottom.  This  can  be  seen  if  the 
normal  mode  is  represented  by  upward-going  and  downward- 
going rays  in  the  layer.  As  the  frequency  of  the  signal  is 
decreased,  the  incidence  ^gracing)  angle  between  the  ray 
and  the  bottom  increase.  Energy  continues  to  be  perfectly 
reflected  and  remains  trapped  in  the  layer  with  no  transfer 
into  the  bottom,  until  the  incidence  angle  becomes  equal  to 
the  critical  angle.  At  that  frequency  the  mode  stops  propa 
gating  in  the  layer  and  energy  is  transferred  into  the 
bottom.  This  concept  can  qualitatively  be  applied  to  a 
tapered  layer  over  a fast  bottom.  As  a propagating  mode 
travels  up  the  wedge,  the  incidence  angle  of  the  equivalent 
ray  increases  with  each  bounce  until  it  reaches  the  oritioa 
angle,  at  which  point  energy  will  be  transmitted  into  the 


bottom. 


Theoretical  and  experimental  studies  of  this  phenomenon 
have  been  carried  out  at  the  NPS  for  three  years  encom- 
passing three  generations  of  experimental  setups. 

Lt.  Edwards  [Ref.  1]  with  his  thesis  adviser.  Professor 
A . B . Coppens  carried  out  initial  investigations  in  a large 
tank  filled  with  brine  solution  over  which  a thin  layer  of 
olive  oil  was  carefully  added.  Since  the  interface  between 
the  two  fluids  was  horizontal,  a wedge-shaped  upper  layer 
was  created  by  inserting  a board  of  pressure- re  lease  material 
at  an  angle  to  the  horizontal.  The  sound  field  was  observed 
in  the  bottom  fluid  at  various  distances  and  depths.  When 
plotted  on  properly  chosen  coordinates , the  data  were 
observed  to  cluster  about  a single  curve  consistent  with 
the  existence  of  narrow  beam  of  sound  enetermg  the  bottom 
at  a relatively  shallow  inclination. 

In  the  meantime.  Professor  Coppens  developed  a simple 
theoretical  model  that  predicted  the  shape  of  the  beam 
pattern. 

In  an  attempt  to  obtain  more  precise  data,  Lt . Netzorg 
l Re f . -),  with  his  thesis  adviser,  Professor  J.V.  Sanders, 
redesigned  the  apparatus  to  consist  of  a wedge  of  water 
separated  from  a small  tank  of  brine  by  a 1-mil  Mylar 
diaphragm  stretched  tightly  within  a frame.  The  bottom  and 
sides  of  the  tank  were  lined  with  absorbing  material  to 
reduce  reflections.  The  results  showed  that  the  entry  point 
of  the  beam  was  close  to  where  it  was  predicted  by  the 
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simple  model.  Agreement  between  experiment  and  theory  was 
fair  but  not  good  enough  to  verify  the  model. 

Subsequent  to  these  first  experiments  a more  refined 
model  has  been  constructed,  with  silicon  oil  in  the  wedge 
and  a large  water  tank  for  the  bottom.  But  the  results 
obtained  in  this  apparatus  by  Professor  Sanders  were  still 
significantly  disturbed  by  what  appeared  to  be  some  sort 
of  interference. 

Experimental  results  from  the  three  generations  are 
compared  to  the  predictions  of  the  simple  model  in  Table  1-1 
and  Fig.  1-1.  Significant  differences  between  theory  and 
experiment  were  observed  on  the  lower  side  of  the  beam 
pattern.  In  one  extreme  example,  two  peaks  of  beam  were 
observed  in  a third  generation  of  experiments.  Generally, 
the  observed  angle  of  depression  of  the  beam  is  shallower 
and  the  beam  width  narrower  than  predioted  by  the  simple 
model . 

At  this  point,  the  research  »ffort  was  redirected  to 
go  back  to  examining  the  amplitude  and  phase  distribution 
along  the  bottom;  these  quantities  provide  the  link  between 
the  acoustic  field  in  the  layer  and  the  formation  of  the 
beam  in  the  bottom. 

The  objective  of  this  research  reported  herein  are  to: 

1)  Develop  a simple,  closed  form  expression  for  the 
pressure  distribution  along  the  bottom  that  can  be 
used  to  reliably  predict  the  beam  pattern  to  the 
bottom. 
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2)  Calculate  the  pressure  (including  phase)  along  the 
bottom  predicted  by  application  of  the  method  of 
images  and  compare  with  experimental  results. 

In  this  report,  two  different  coordinate  systems  are  used. 

The  first  is  a special  coordinate  system  applicable  to  the 
simple  model;  the  origin  is  located  at  a point  on  the  bottom 
where  the  lowest  mode  reaches  cutoff.  The  second  is  a 
general  coord.nate  system  more  suitable  for  the  image 
calculations  and  experiment;  the  origin  is  located  at  the 
apex.  These  two  coordinate  systems  are  shown  in  Fig.  I-2a,b. 


>mpa.rl son  of  the  experiments  with  the 
simple  model 


II.  THEORY 


A.  PRESSURE  DISTRIBUTION  PREDICTED  BY  A SIMPLE  MODEL 
1.  Basic  Concept 

The  purpose  of  this  part  is  to  summarize  the 
theoretical  background  necessary  to  conceptualize  the 
propagation  of  sound  in  a tapered  fluid  layer  above  a fast 
bottom,  using  a combination  of  normal  modes  and  ray  theory, 
a.  Layer  of  Uniform  Depth 

One-dimensional  propagation  in  a fluid  layer 
of  uniform  depth  H with  speed  of  sound  c^  and  density  p^, 
overlying  a second  layer  with  a faster  speed  of  sound 
c^,  and  density  p^  can  be  exactly  represented  by  an  equation 
of  the  form 

i (u)t-Kx,x'  ) 

P > [ A sin ( K , z')e  n (II-l) 

n 

n 

where  z'  is  the  depth,  x'  is  the  horizontal  distance, 

K , and  K , are  the  components  of  the  propagation  vector 

2 „ X 

n n 

k for  the  n*"*1  mode  in  the  z'  and  x'  direction  respectively. 
An  is  the  magnitude  of  the  nl  mode.  For  propagating  modes, 
the  propagation  vector  Tc  can  be  represented  as  in  Fig.  II-l. 

— t 

The  angle  between  k and  the  horizontal  is  ■? . Snell's  law 
gives  the  relationship  between  the  speeds  of  sound  in  the 
layer  and  the  critical  angle  $c . 
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Fig.  II-l.  The  propagation  vector  k and  its  components 


0 

c 


arc  cos 


(II-2) 


As  long  as  0 < 9c  the  energy  is  totally  reflected.  How- 
ever, when  9 > 9C  energy  is  transmitted  through  the 
previously  totally-reflecting  surface. 

To  satisfy  the  requirements  of  propagation, 

k , must  be  a constant,  real  number  for  each  mode  and  from 
n 

Ref.  3, 


where 


tan  k , H 
z'n 


(II-3) 


( II-4 ) 


At  cutoff  0 = 9c  and  then  the  geometry  of  Fig.  II-l 
implies  k , = — so  that 

x c 2 
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(n  - j)  it 


(II-5) 


k , H 
2 n 


and  for  a mode  with  frequency  just  above  its  cutoff 
frequency  f there  is  a pressure  antinode  at  the  bottom. 
Now  by  writing  Eq.  II-l  in  the  form 


. ik  , z' 
A z 

r (e  n 

1 2l  (e 

n 


a lkz'nZ  j ei (wt-kx' ) 


i l i An<, 

n 


i(“)t-k  , x'+k  , z') 
n n 


- e 


i (uit-k  , x'  -k  , z ' ) 
x n 2 


n } 


(II-6) 


we  have  two  plane  waves  that  propagate  with  wave  number 

^ A A 

k * k «:t  ± k iz'.  If  we  assume  that  the  frequency  is 

n xn 

slowly  lowered,  higher  modes  are  cutoff  as  the  frequency 

I TT 

decreases  with  k , approaching  (n  - y) — 

z n i « 

b.  A Tapered  Layer 

Let  us  now  consider  a wedge  with  angle  6 as 
indicated  in  Fig.  II-2  and  II-3,  where  is  the  distance 


0 to  the  apex,  and  x'  is  the  distance  of  closest 
c s 


from  6 

approach  (Turnaround  distance)  of  the  indicated  ray.  In 
this  case  each  bounce  off  the  bottom  increases  0 by  23.  At 
the  point  where  0 = 0c  energy  begins  to  penetrate  into  the 
lower  layer.  The  resulting  attenuated  ray  bounces  off  the 
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top  and  makes  another  collision  with  the  bottom,  losing 

still  more  energy.  This  process  continues  until  the  wave 

reaches  the  "Turnaround  distance"  x'  defined  by: 

s J 


X ' » X'  (1  - sin  6 ) 

so  c 


(II-7) 


For  the  wedge  then,  the  frequency  remains  constant  while 

d increases  with  each  bounce  off  the  bottom  until  that  depth 

H is  reached  for  which  6=9  where 
O c 


(n  - -j)  tt 


From  Fig.  II-l  we  obtain  at  cutoff 


( n-8 : 


= k sin  9 


( II-9 ) 


where 


(ii-io: 


and  substituting  Eq.  (II-9)  into  Eq.  (II-8)  yields 


(n  - |)  a 


(11-11) 


Knowing  H allows  calculation  of  the  value  of  x'  = X'  where 
° o 

energy  will  begin  to  be  transmitted  into  the  bottom. 
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!•  Pressure  Distribution  along  the  Bottom 

A.B.  Coppens  [Ref.  4]  formulated  an  approximate 
model  to  predict  the  characteristics  of  the  transmitted 
beam  in  the  bottom.  The  expressions  for  the  pressure  ampli- 
tude and  are: 

, „ -A 

A ■ AQ  e (11-12) 


/- *r~  / sin  9 7 

*«■  ■ ■ t-J1  - ‘nf1 

o 

where  x' , X^,  are  defined  from  the  geometry  of  Fig.  II-l, 
^rora  E(3*  Ac  is  the  pressure  amplitude  at  x'  =0, 


c 


( 11-16 ) 


i 


- cos  6 

c 


fl 

c2 


di-17) 


x 1 1 2 

In  the  case  where  jrr  j cos  8^  this  model  may  be  simpli- 

o 

fied  further 


* 


°1  tan,,c  2.?,x',3/2 
2 o 


,11-18) 


and 


a 


sm 
1 - 


(11-19) 


e 


arc  sin 


sin 
1 - 


(11-20) 


3 . Pressure  Distribution  under  the  Bottom 

In  this  part  we  present  without  detailed  derivation 
A.B.  Coppens  [Ref.  21  work  concerning  the  pressure  distri- 
bution under  the  bottom  based  on  the  simple  model  and 
obtained  through  the  use  of  the  Green's  function  and  the 
method  of  stationary  phase.  Two  basic  assumptions  were 
considered  for  this  development. 

a.  The  plane  wave  has  wavefronts  that  are  parallel 
to  the  shore  line. 


H 


b.  The  pressure  distribution  on  the  bottom  can  be 
viewed  as  an  acoustical  source  which  generates 
the  beam  of  sound  projected  into  the  bottom. 

The  geometry  is  indicated  in  Fig.  II-4. 

The  Green's  function  formulation  for  this  case  can  be 
developed  through  the  geometry  of  Fig.  II-4  as  follows: 

At  the  point  (x',z‘)  the  instantaneous  pressure  is  given  by 

P = P(r',t)  = P(r')e1L)t  (11-21) 


where 


f'r'>  - I ?s(x">!f’-|2’=o  dx"  (II-22) 

—co 

The  integrand  of  Eq.  (11-22)  is  formed  from  the  pressure 
distribution  at  the  bottom 


Pg(x") 


/ sin©  0 

_i  t/1'(r^r)  x" 

A(x")  e Xy  L'X^ 


1 - 


(II-23; 


and  the  Green's  function  for  a line  source  is 


G(r*  ,r")  = 


i tt  H (2)  ( |r'-r"|  ) 
o 


(11-24) 


23 


1 


In  Eq.  (11-23)  the  term 


is  an  a postiori  correction 


1 - 


x 

X1 


' o 

for  the  cylindrical  convergence  of  sound.  In  the  absence 
of  losses  into  the  bottom  and  if  there  is  no  energy  loss 
within  for  the  wedge,  the  pressure  amplitude  would  behave 
as 


(11-25) 


also  from  the  geometry  of  Fig.  II-2  we  have 


o 


so  that 


(11-27) 


The  term  A(x")  is  the  pressure  amplitude  on  the  bottom 
and  the  exponential  term  is  the  corresponding  phase. 

This  pressure  distribution  is  an  integral  of  the 
general  form 


f F(x)  eG[K)  dx  (11-28) 


• 00 


25 
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where  F(x)  varies  slowly  compared  with  the  exponential 

factor.  The  value  of  the  integral  can  be  well  approached 

with  the  method  of  stationary  phase. 

Using  the  above  method  and  assuming  the  simple 
x"  1 2 

model  for  <<  cos  6c>  and  for  small S the  following 
o 

expression  for  the  pressure  distribution  was  developed: 


Pi  ^c  _1_  3 

P2  S 6/2 


(11-29) 


where 


/ 2(1  - cos  9 ) _ sin  9 

tan  6 ~ tan  9, 


(11-30) 


2 . Pressure  Distribution  along  the  Bottom  Predicted 
by  the  Method  of  images 

Since  the  simple  model  did  not  yield  results  agreeing 
with  the  experimental  measurements,  it  was  decided  to  use 
the  exact  but  more  cumbersome  method  of  images  to  attempt 
to  predict  the  pressure  distribution  on  the  bottom. 

Let  the  source  be  a distance  X from  the  apex  of 
the  wedge  and  a depth  H below  the  sea  surface  as  indicated 
in  Fig.  II-5. 

Define 


A . H 

|>d  = arc  sin  ^ 


(II-31) 


L 


i 


1 


— 


and 


* 


o 


4 


(11-32) 


where  8 is  the  angle  between  the  sea  surface  and  bottom, 

4>d  is  th®  angle  formed  at  the  apex  between  the  source  and 
the  sea  surface,  and  is  the  analogous  angle  for  the  source 
and  the  bottom. 

If  4>n  is  the  angle  formed  at  the  apex  between  the 
n image  and  the  bottom  then  Fig.  I 1-5  reveals 


= tJ  + * 26  - 4> 


= + 2(6“>t>  ) = 2S+it> 

zoo  o 


4>2  * 34)^  + 4(6->J>0)  = 46-4>0 


( 11-33) 


or  more  generally 


* (n+l)4>o  + n(6“4i0)  = n6  + $Q  for  n even 

* nit  + (n+1)  (6~4>  ) = (n+l)6  - 4>  for  n odd 

no  o o 

(11-34) 
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This  set  of  equations  can  be  reduced  further  to  one 
expression : 

*n  * tn  * |[1-  (-1)  n]  }g  + (-l)n4>0  - 2 INT  l^)  8+  (-1 ) %0 

( II-35a) 

Also  for  the  symmetric  member  n'  of  the  paired  images 
we  obtain 


{ n 1 + ill-  (-1) 


} (3  - ( - 1 : 


2 INT[2-t^]i*  +(-l)n' 


(II-35b) 


where  INT ( 1 denotes  the  largest  integer  which  is  equal 

to,  or  smaller  than  the  argument. 

Considering  the  geometry  of  Fig.  II-5  and  II-6  we 

can  derive  expressions  for  the  distances  r and  r , which 

n n 

represents  the  distance  between  a receiving  point  on  the 
bottom  at  distance  x from  the  apex  and  the  nfc^  and 
image  respectively.  Also  we  can  derive  expressions  for 
the  angles  f'nQ  and  8 ,q  which  formed  at  the  receiving  point 
between  the  bottom  and  the  images  n and  n'  respectively. 
These  expressions  are: 
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Fig. 


II-6. 


Geometric  development  of  the  image's  position 
with  respect  to  the  receiving  point 
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arc  tan  [ 


(II-37a) 


sin 


n 


no 


COS  0 - - 

n x 


sin 


0 , * arc  tan  [ 

no 


^n ' 


cos 


n ' 


xJ 

X 


(II-37b) 


Let  us  now  define  9 and  9 , as  the  angles  of 

nm  n m r 

incidence  for  the  m bounce  from  the  bottom  for  the  n 

and  n'  image  respectively,  m = 1,2,3,  ...  (The  oth  bounce 

is  the  last-one  before  reaching  the  receiver.)  From  the 

geometry  of  Fig.  II-5  and  Fig.  II-7  we  obtain  expressions 

for  9 as  follows: 
nm 
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The  general  expression  is 


0 = 0 - 2mS 

nm  no 
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Also  for  the  symmetric  member  n'  of  the  nfc^  pair  we  obtain 


0 * 0 , - 2m6 

n m no 


1,1.  i,  u 

and  the  maximum  number  of  bounces  of  the  n and  n' 


( II-39b) 


image  is 
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If  the  receiving  point  is  on  the  bottom,  then  from  Fig.  II-S 


we  obtain: 


1.  The  distance  from  each  member  of  a pair  of 
images  (n,n')  to  the  receiving  point  is  the 
same,  i.e. , r * rn • • 

2.  The  number  of  reflections  ;from  the  bottom) 

t h 

from  the  n'  image  to  the  receiving  point  is 

th 

one  more  than  the  n image. 

3.  The  grazing  angle  to  the  bottom  is  the  same 
for  each  corresponding  bottom  bounce  so  that 
the  reflection  coefficients  for  corresponding 
bounces  are  equal,  except  the  last  bounce  for 
the  n'  ‘ imaue  . 


Consequently  Rnm  - ?n.m  where  and  Rn>m  are  the  reflec 

tions  coefficients  for  the  bounce  of  the  n and  n' 

image  respectively. 

t hi 

The  reflection  coefficients  for  the  n image  are 


Let  us  now  consider  the  pressure  P(x)  at  a point  x on  the 
bottom.  The  contributions  from  the  paired  images  are 
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and 


N , -ikr  , M INT  [2^1] 

i lb e n !rf-  <-l> 
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n=0 


m=0 


(The  nfck  and  n,fc^  images  differ  only  in  that  nfc^  image  has 

no  R „ reflection.) 

—no 

The  complex  pressure  is 


N -ikr  INT [— y— ] , M 

P(X)  * l {-=-  e n (-1)  (l+p^-)  n Rnm> 

- r R _ -nm 

n—0  n -no  m—  0 (11  — 47) 


A computer  program  was  constructed  to  calculate  this 
complex  pressure  (amplitude  and  phase) . This  program, 
with  the  corresponding  flow  chart,  is  presented  in  Appendix  A. 

In  case  of  a source  so  distant  that  it  can  be 
assumed  to  be  at  r * * then  the  geometry  of  Fig.  II-9 
yields  the  approximation 
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0 ~ 4) 

no  *n 


(11-48) 


which  simplify  the  mathematics  significantly.  Here 


?no  * ?(eno)  = ?(V 


(11-49) 


In  this  case,  the  expression  for  reflection  coefficients 
are  (see  Appendix  B) 
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for  n > 2 


where  the  right  hand  sides  of  Eqs.  (11-50)  depend 
only  on  n,  then  define 
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(11-51) 

r(n-2)Rn 

for 

n > 

2 

Also,  the  geometry  of  Fig.  II-9  gives 

r = r - x cos  ? (11-52) 

n n 

Substitution  of  Eq.  ((11-50),  (11-51),  (11-52)  into 

Eq.  (11-47)  yields 
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The  last  expression  for  pressure  is  the  approximation  for 
a distant  source,  which  is  described  in  detail  elsewhere 
[Ref.  4). 

The  advantage  of  Eq.  (11-53)  with  respect  to  Eq.  (11-47) 
is  that  because  the  elimination  of  product  term,  the 
calculations  are  greatly  simplified. 
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III.  EXPERIMENT 


Initially  it  was  intended  that  the  beam  pattern  in  the 
bottom  would  be  exhaustively  investigated,  and  a larger 
experimental  system  with  silicon  oil  in  the  wedge  and  a 
large  water  tank  for  the  bottom  was  designed.  But,  even 
in  the  new  system,  it  was  observed  that  the  beam  pattern 
was  being  significantly  disturbed  by  unknown  interferences 
that  caused  the  apparent  angle  of  the  beam  to  move  to 
shallower  angles  than  predicted  by  the  simple  theory  des- 
cribed in  II-A  and  made  the  beam  seem  narrower  than  predicted. 
This  interference  was  definitely  not  from  the  walls  or 
bottom  of  the  water  tank.  It  might  have  been  associated 
with  the  wedge,  but  whether  it  was  an  inherent  property  of 
wedge  propagation  omitted  from  the  theory  II-A  or  merely 
an  artifact  of  the  experimental  setup  was  the  question. 

In  an  attempt  to  study  this  question  more  thoroughly, 
the  experiment  was  redirected  to  investigate  the  amplitude 
and  phase  of  pressure  on  the  bottom  of  the  wedge. 

A.  APPARATUS 

The  experimental  facility  consisted  of  a wedge  of  sili- 
con oil  overlying  a large  tank  of  fresh  water;  the  two 
fluids  being  kept  apart  by  a thin  Mylar  membrane.  A photo- 
graph of  the  apparatus  is  shown  in  Fig.  III-l.  The  sound 
source,  positioned  at  the  end  of  the  silicon-oil  wedge, 
sent  a pulsed  signal  toward  the  apex  of  the  wedge.  It  was 
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possible  to  position  the  receiver  anywhere  within  the  fresh 
water  region,  but  in  this  experiment  the  receiver  was  mainly 
placed  clos*.  to  the  bottom  of  the  wedge.  Tne  geometry  is 
shown  in  Fig.  III-2. 

1.  Water  Tank  (Fig.  III-3a,b,c) 

The  water  tank,  constructed  of  3/4  in.  thick 
plywood,  had  the  inside  surfaces  painted  with  Varathane 
but  were  otherwise  uncoated.  It  was  filled  with  tap  water 
at  room  temperature. 

2.  Wedge  (Fig.  III-4) 

The  silicon  oil  was  contained  in  a frame  whose 
sides  were  made  of  1/2  in.  x 2 in.  aluminum  and  whose 
bottom  consisted  of  a thin  Mylar  sheet.  Half-mil  thick, 
aluminum-coated  Mylar  was  fastened  to  the  bottom  of  the 
frame  by  first  taping  a sheet  of  Mylar,  metalized  side  up 
to  a flat  surface.  Then  the  aluminum  frame  was  placed  on 
the  Mylar,  and  the  edges  outlined  by  a sign  pen.  A bead 
of  "Royal  bond  GRIP"  contact  cement  was  applied  to  the 
bottom  of  the  frame  and  to  the  marked  areas  of  Mylar.  After 
the  cement  was  completely  dry  (about  15  minutes) , the 
frame  was  placed  on  the  Mylar  and  about  50  pounds  of  weight 
placed  on  top  of  the  frame.  After  the  contact  cement  had 
cured  for  12  hours,  a coat  of  "Silver  print  conductive  paint" 
was  applied  to  the  outside  joint  between  the  Mylar  and  the 
frame.  This  provides  electrical  contact  between  the  Mylar 
and  the  frame.  (This  was  necessary  to  allow  measurement  of 
the  slope  of  the  wedge  as  will  be  explained  in  a later 
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i 

section.)  The  excess  Mylar  was  trimmed  flush  with  the 
outside  edge  of  the  frame  by  a razor  blade.  As  a final 
step,  a thin  coat  of  the  "Silicone  rubber  sealant"  was 
applied  to  both  the  outside  and  inside  of  the  Mylar-frame 
interface  to  prevent  leaks.  Because  the  tension  in  the 
Mylar  tended  to  bend  the  walls  of  the  frame  inward,  a 
3/8  in.  diameter  aluminum  rod  was  used  as  a "frame  spreader" 
to  spread  the  frame  and  keep  Mylar  flat.  The  length  of 
the  frame  spreader  was  adjustable.  It  was  placed  across 
the  shorter  dimension  of  the  frame  and  adjusted  until  the 
Mylar  film  was  mirror  flat.  Aluminum  arms  at  each  end  were 
used  to  suspend  the  wedge  from  the  top  of  the  water  tank. 

The  longer  arm  was  connected  to  the  apex  side  of  the  wedge 
and  extended  in  the  propagated  direction.  The  slope  of  the 
wedge  was  set  by  raising  the  end  of  this  arm  with  wooden 
blocks.  The  other  arm,  perpendicular  to  the  propagation 
direction,  was  connected  to  the  source  side  of  the  wedge. 

The  horizontal  level  of  the  wedge  was  adjusted  by  means  of 
screws  at  the  ends  of  this  arm. 

3 . Source  Transducer 

The  source  transducer  was  made  by  Lt.  Netzorg. 

Construction  details  are  discussed  in  his  thesis  [Ref.  2J. 

The  design  of  the  transducer  was  quite  well  done:  it  has 
sufficient  horizontal  directivity  so  that  sound  is  not 
reflected  from  the  side  of  the  frame,  there  is  negligible 
sound  propagation  in  the  backward  direction,  and  its  active 
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face,  2 cm  high  by  8 cm  wide,  is  the  largest  size  allowable 
in  the  wedge.  The  experimentally  measured  directivity  is 
very  close  to  the  directivity  predicted  for  a rectangular 
source  (Fig.  III-5a,b,c) . The  theoretical  directivity 
is  described  as  follows: 


H(0,4>) 


sin (^kL^sinO ) 
•^kL^sine 


sin  ( jkL^sin«p ) 
ikL0sin>J> 


(Ref.  3,5] 


(III-l) 


where 

= Horizontal  dimension  of  an  active  face 
of  a rectangular  source. 

L-,  = Vertical  dimension  of  an  active  face  of 

a rectangular  source. 

K = £ = rJL  wave  number 

c \ 

0 = Horizontal  angle  measured  from  a major 

radiation  axis. 

-S  = Vertical  angle  measured  from  a major 
radiation  axis. 

As  a matter  of  convenience,  the  directivity  measurements 
were  done  in  the  water  instead  of  silicon  oil.  The  fre- 
quency used  in  these  measurements  was  adjusted  to  give  the 
same  wavelength  as  produced  when  the  source  operates  at  its 
prescribed  frequency  in  silicon  oil.  A frequency  of 
127.8  kHz  in  fresh  water  ( c - 1503.8  m/sec)  produces  the 
same  wavelength  as  100  kHz  in  a silicon  oil  (c^  = 1176.5  m,  sec) . 
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For  more  accurate  comparison  between  the  experimental 
results  and  the  method  of  images,  the  vertical  directivity 
was  employed  in  the  computer  model  "WEDGE  0"  (APPENDIX-A) . 


The  source  transducer  was  mounted  on  a support  that 
allowed  it  to  be  located  anywhere  in  the  upper  fluid  layer. 
Vertical  adjustments  were  made  by  means  of  a micrometer  with 
1/2  in.  of  travel. 

4.  Receiver  and  Attachment  (Fig.  III-6a,b) 

A LC-10  hydrophone,  manufactured  by  Celesco  Indus- 
tries, was  used  as  the  receiver.  It  was  possible  to  locate 
the  receiver  anywhere  in  the  fresh  water.  Since  the  main 
purpose  of  the  experiment  was  to  investigate  the  sound 
pressure  distribution  on  the  bottom,  the  hydrophone  was 
usually  located  as  close  to  the  bottom  as  possible.  The 
attachment  mounted  on  the  frame  of  the  wedge  as  shown  in 
Fig.  III-6a. 

At  the  beginning  of  this  thesis,  some  measurements 
were  made  of  the  beam  pattern  in  the  lower  fluid.  Then  a 
different  attachment  was  used  (Fig.  III-6b)  which  was  mounted 
on  the  large  arm. 

In  both  cases,  the  depth  of  they  hydrophone  was 
adjustable.  In  all  cases,  measurements  were  made  in  the 
vertical  plane  that  contained  the  acoustic  axis  of  the  source 
transducer. 

5.  Slope-measurement  Assembly  (Fig.  III-7a,b,c) 

The  slope  measurement  assembly  consisted  of  a pin, 
micrometer,  and  attachment.  Fine  vertical  adjustments  of 


44 


the  pin  were  made  by  means  of  the  micrometer  with  1/2-in. 
of  travel.  The  attachment  was  mounted  on  the  wedge  frame 
and  could  be  located  anywhere  in  the  wedge.  The  pin  was 
electrically  insulated  from  the  frame,  so  that  when  an 
Ohmmeter  was  attached  between  the  frame  and  the  pin,  the 
precise  depth  at  which  the  pin  contacted  the  metalized 
Mylar  could  be  determined. 

6 . Electrical  Setup 

Figure  III-8  is  a block  diagram  of  the  electrical 
setup  which  consists  of  a driving  system,  a receiving  system, 
and  a measuring  system.  The  transducer  was  driven  by  a 
General  Radio  Oscillator  (Type  1310)  and  Tone  Burst  Generator 
(Type  1396-A) , polarized  by  a 300  Volts  D.C.  biasing-voltage 
[Ref.  2] . The  driving  frequency  (100  kHz)  was  measured  with 
a Hewlett  Packard  electronic  frequency  counter  (Type  5233L) 
and  the  driving  waveform  displayed  on  one  channel  of  a 
Tektronic  Dual-Beam  Oscilloscope  (Type  565) . The  received 
signal  was  amplified  (40  dB)  by  a Hewlett  Packard  Amplifier 
(Type  465-A) , and  band  passed  (40  kHz  ~ 200  kHz)  through  a 
Spencer-Kennedy  Laboratories  Band  Pass  Filter  (Serial  1683 
Model  302) . The  received  signal  was  displayed  on  the  other 
channel  of  the  Tektronix  Dual-Beam  Oscilloscope  (Type  565) 
and  the  amplitude  was  measured  on  this  display.  The  received 
signal  also  went  to  a Hewlett  Packard  Oscilloscope  (Model 
120B)  where  its  phase  was  compared  to  that  of  the  driving 
signal . 


7 . Auxiliary  Apparatus 

To  measure  the  speeds  of  sound  and  densities  of 
fluids,  other  pieces  of  equipment  were  used. 

Transducer 

Two  Hydrophones  LC5-2 

Analytical  Balance 

ASTM  Vessel 

B.  PROCEDURE 

Before  beginning  the  experiments,  the  speeds  of  sound 
and  densities  of  the  fluids  (a  silicon  oil  for  the  upper 
layer,  fresh  water  for  the  bottom)  and  the  slope  of  the 
bottom  was  measured  with  high  accuracy. 

1.  Sound  Speed  Measurements  (Fig.  III-9) 

The  speeds  of  sound  were  measured  in  a pan  which 
had  dimensions  10  3/4  in.  x 7.0  in.  x 1 3/8  in.  About  1500  cc 
of  a test  fluid  was  carefully  poured  into  a pan  to  prevent 
the  entrapment  of  air  bubbles.  Two  LC5-2  receivers  and  a 
transducer  were  used  to  measure  the  speeds  of  sound.  The 
distance  between  two  receivers  was  set  at  exactly  20.0  cm. 

And  the  transducer  was  driven  by  a 150  kHz  pulsed  signal. 

The  time  of  transit  for  one  receiver  to  the  other  was  measured 
using  the  time-delay  feature  of  the  dual-beam  oscil loscope . 

The  speed  of  sound  was  then  calculated  by  dividing  the 
distance  by  the  transit  time. 

2 . Density  Measurments 

The  densities  of  the  two  fluids  were  found  by  first 
weighing  a 100  ml,  ASTM  vessel  on  an  analytical  balance 
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readable  to  one-tenth  mg.  The  vessel  was  then  filled  to 
the  100  ml  mark  and  weighed  again.  By  subtracting  the  first 
weight  from  the  second  weight,  the  weight  of  100  ml  of  the 
test  fluid  was  determined.  The  weight  of  the  test  fluid  in 
grams  divided  by  the  100  ml  volume  and,  the  quotient  times 
1000,  gave  the  density  of  the  test  fluid  in  kg/m-*. 

3 . Measurement  of  the  Slope  Angle 

The  approximate  slope  was  set  by  raising  the  end 
of  the  long  arm  a precalculated  height  with  the  use  of 


wooden  blocks.  Then  3500  cc  of  silicon  oil  were  slowly 
poured  into  the  wedge,  and  the  frame-spreader  was  adjusted 
to  keep  the  bottom  flat.  The  water  level  in  the  tank  was 
then  adjusted  to  keep  the  apex  line  (shore  line)  straight. 

The  source  transducer  was  set  at  about  middle  depth  on  the 
wedge  centerline.  The  slope  angle  was  accurately  measured 
with  the  slope-measurement  assembly  described  in  the  appara- 
tus section.  Because  the  top  and  bottom  of  the  frame  of 
the  wedge  are  parallel,  if  the  Mylar  sheet  is  flat  the 
distance  from  a reference  level  to  the  bottom  should  be  the 
same  at  all  points.  The  distance  from  a reference  level  to 
the  surface  and  to  the  bottom  were  measured  at  several  points 
on  the  wedge  centerline,  and  the  data  were  plotted.  The 
slope  was  calculated  from  the  plot  as  shown  in  the  example 
illustration  in  Fig.  III-10. 

4 . Measurement  of  Pressure  Amplitude  along  the  Bottom 
The  electrical  setup  of  the  measurements  was  shown 
in  Fig.  III-3.  The  oscillator  output,  set  at  100  kH2,  was 
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fed  into  the  tone  burst  generator.  A pulsed  signal  of 
36  cycles  of  the  100  kHz  input  with  a closed-qate  duration 
time  of  10  msec  (repetition  rate  was  about  100  Hz)  went 
through  the  transducer  biasing  network  and  then  to  the 
source  transducer.  The  driving  signal  was  displayed  on  one 
channel  of  the  dual-beam  oscilloscope  and  adjusted  for  a 
peak-to-peak  voltage  of  10  Volts.  The  receiving  signal  was 
displayed  on  the  other  channel  of  the  same  dual-beam  oscillo- 
scope. This  was  very  convenient  because  the  receiving  sig- 
nal was  always  compared  with  the  driving  signal  on  the  same 
oscilloscope.  Reflections  of  the  receiving  signal  were  not 
observed. 

5 . Measurement  of  Pressure  Phase  along  the  Bottom 

The  phase  of  the  receiving  pressure  was  compared  to 
that  of  the  driving  signal  on  an  oscilloscope . Fig.  III-8 
also  shows  the  electrical  setup  for  these  measurements.  The 
signal  from  the  oscillator  was  directly  fed  into  the  x-axis 
of  the  oscilloscope.  The  receiving  signal  was  amplified 
and  filtered,  then  fed  into  the  y-axis  of  the  same  oscillo- 
scope. The  signal  on  the  x-axis  was  continuous  with  a con- 
stant amplitude  (adjusted  to  10  Volts  peak-to-peak  voltage) . 
On  the  other  hand,  the  signal  applied  to  the  y-axis  was 
pulsed  and  had  an  amplitude  that  depended  on  the  receiver 
position.  The  pattern  observed  on  the  oscilloscope  was  an 
ellipse  crossed  by  a line  on  the  x-axis  (like  "0"  of  the 
Greek  alphabet) . The  line  occurred  during  the  time  the 
received  pressure  was  zero,  i.e.,  between  pulses.  When  the 


signal,  generated  by  the  oscillator,  was  sent  directly  to 
the  transducer,  bypassing  the  tone  burst  generator,  the 
pattern  on  the  oscilloscope  was  the  same  ellipse  but  without 
the  line  on  the  x-axis.  This  indicates  that  the  pulse  is 
sufficiently  long  to  allow  accurate  phase  determinations. 

(It  also  showed  that  reflections  in  the  tank  were  negligible.) 
But  measurements  were  still  difficult,  because  the  amplitude 
of  the  receiving  signal  was  not  constant  and  the  phase 
changed  very  quickly  with  changes  in  the  receiver  position. 

The  phase  measurements  were  done  independent  of  the  ampli- 
tude measurements.  The  position  of  the  receiver  was  recorded 
when  the  phase  was  0°,  90°,  180°  and  270°.  (Corresponding 
to  a straight  line  with  positive  slope,  a symmetric  ellipse 
and  a straight  line  of  negative  slope.) 
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Fig. ///-7b  Slope-measurement  assembly  and  wedge 
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b'igJII-9  Liound  speed  measurements 


DISTANCE  FROM  THE  SOURCE  (cm) 
Fig.///-IO  Measurement  of  the  slope  (CASE  2) 


IV.  RESULTS  AND  COMPARISON 


Measurements  were  conducted  under  two  different 
conditions : 

Common  Parameters 

Speed  of  sound  in  the  wedge  c.  = 1176.5  m/sec 

(silicon)  i 


Speed  of  sound  under  the 
(fresh  water) 

Density  of  a silicon  oil 

Density  of  a fresh  water 

Frequency 

Pulse  length 

Closed  gate  duration 

Pulse  repetition  rate 

CASE  1 

Slope  angle 

Source  distance  from  the 
Source  depth 


bottom  c^  = 1503.8  m/sec 

PL  = 1015.0  kg/m3 

= 995.0  kg/m3 

f = 100  kHz 

Lp  = 36  wavelengths 

t =10  msec 
c 

f = 100  Hz 
P 

6 =0.047  rad 

apex  X = 7 0 . 0 cm 

H = 1 . 6 4 cm 


CASE  2 

Slope  angle 

Source  distance  from  the  apex 
Source  depth 


8 = 0.042  rad 

X = 7 5 . 5 cm 
H = 1 . 5 7 cm 


The  speeds  of  sound  and  the  densities  of  the  fluids  were 
measured  at  18.5°C.  The  two  cases  differed  mainly  in  the 
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slope  angle:  for  the  first  case  3 = 2.7°  and  for  the 
second  S = 2.4°.  These  are  fairly  realistic  values: 
according  to  Sverdrup- Johnson' s - Fleming  [Ref.  6],  con- 
tinental slopes  off  mountainous  coasts  have,  on  the 
average,  a slope  of  about  3.5°,  whereas  off  coasts  with 
wide,  well-drained  coastal  plains,  the  slope  is  about  2°. 
However,  there  are  many  different  estimates  based  on  differ- 
ent measurements  at  different  locations.  For  extreme 
examples,  estimates  of  the  average  continental  slopes  range 
from  about  5°  to  about  0.1°  [Ref.  7], 

The  distance  from  the  apex  to  the  source  was  limited 
to  less  than  about  6 times  the  lowest  possible  cutoff 
distance  xQ  (xq  = 9.92  cm  for  CASE  1,  xQ  = 11.34  for  CASE  2). 

A.  EXPERIMENT  VS.  SIMPLE  MODEL 
1.  Pressure 

Figures  IV-la,b  show  the  experimentally-obtained 

pressure  distribution  along  the  bottom.  In  both  cases 

there  was  a dominant  pressure  peak  near  the  lowest  cutoff 

point  but  shifted  slightly  in  the  direction  away  from  the 

apex.  To  compare  the  simple  model  with  the  experiment, 

data  and  the  theoretical  curve  were  plotted  on  the  same 

graph  for  the  interval  0 <_  x <_  2xq  (Fig.  IV-2a,b)  . There 

are  significant  differences  between  the  experimental  results 

and  the  simple  model.  The  theoretical  curves  decay  quickly, 

and  the  model  predicts  that  there  should  be  no  energy 

transmitted  into  the  bottom,  outside  the  region  x < x < X 

s — — o 
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(the  small  amount  observed  in  the  simple  model  results 
from  the  approximations) . On  the  other  hand,  the  experi- 
mental curves  have  a long  tail  with  energy  transmitted 
continuously  into  the  bottom  all  the  way  to  the  apex,  and 
a pressure  peak  shifted  significantly,  so  that  energy  is 
leaked  into  the  bottom  before  the  cutoff  depth,  x > xQ. 

2 . Phase 

The  phase  distribution  was  measured  in  the  region 
0 £ x £ 6xq  for  CASE  2.  The  rate  of  change  in  the  measured 
phase  was  almost  constant  but  with  small  irregular  fluctua- 
tions. Figure  IV-3  compares  the  experimental  results  within 
the  region  0 £ x £ 2xq  with  the  results  of  the  simple  model 
calculated  from  Eg.  (11-13) 


kx' 


Phase  = -k  , x ' = 
x' 


•x„  -r-(i  - £-) 

O C ^ Xo  1 1 


sm  a 

V 


(IV-1) 


For  comparison  with  the  experiment,  the  phase  was  restricted 
to  the  interval  -it  to  +tt. 

Figure  IV- 4 shows  the  apparent  angle  of  incidence 

sin  0 sin  0 

6 = arc  sin  ( -x)  = arc  sin  ( — — ) 

^ X X 

Xo  Xo  (IV -20) 
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I 


I 


L 


which  is  calculated  from  the  phase  (-k^.x')  as  discussed 
in  Section  II-A2.  The  graph  shows  that  comparison  between 
the  simple  model  and  the  experiment  is  better  accomplished 
by  studying  0 . The  rate  of  change  in  the  phase  is  almost 
constant  when  x >_  xQ,  but  the  incident  angle  9 changes 
significantly  in  this  region.  The  incident  angle  6 increases 
with  decreasing  distance  x.  At  x - xq  the  incident  angle 
0 becomes  the  critical  angle  0c  and  Eq.  (IV-1)  shows  that 
-kx,x'  = 0.  Subsequently  the  incident  angle  increases  more 
rapidly  as  x decreases  from  xQ  to  xg.  At  the  turnaround 
point  x = x the  angle  9 becomes  90°  and  again  -k  ,x'  = 0. 

G X 

Thus,  energy  can  be  transmitted  into  the  bottom  only  for 
xg  <_  x <_  xq  according  to  the  model.  Precise  experimental 
investigation  of  the  incident  angle  was  difficult  because 
the  rate  of  change  in  the  phase  was  very  small.  Significant 
differences  between  experiment  and  the  simple  model  were 
observed.  From  the  experiment,  there  is  no  turnaround 
point.  The  experimental  pressure  distribution  has  a long 
tail  which  means  energy  can  transmit  into  the  bottom  all 
the  way  to  the  apex.  And  the  slopes  of  the  phase  curves 
(Fig.  IV-3 ) show  significant  differences  between  experiment 
and  the  simple  model. 

B.  EXPERIMENT  VS.  METHOD  OF  IMAGES 

The  computer  program  "WEDGE  0"  was  designed  for  a com- 
parison between  experiment  and  the  method  of  images.  The 
program  has  the  flexibility  of  an  adjustable  source  distance 
and  inclusion  of  the  vertical  directivity  of  the  source. 
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1.  Pressure 


Before  beginning  the  comparisons,  the  computed 
pressure  distributions  were  observed  with  several  combina- 
tions of  distance  (x  = 100  m)  or  near  source  (x  = 1 m)  and 
with  directional  or  omni-directional  source.  In  all  cases 
the  major  pressure  peaks  and  tails  were  similar;  the 
distance-source  case  showed  slightly  smoother  pressure 
distributions  and  was  little  effected  by  source  directivity. 

On  the  other  hand,  directivity  becomes  important  if  the 
source  was  near.  As  a result,  computations  were  done  with 
the  actual  source  distance  (x  = 70.0  cm  for  CASE  1, 
x = 75.5  cm  for  CASE  2)  and  source  directivity  described 
in  Section  III-A3,  and  these  results  were  compared  with 
the  experiment. 

Figures  IV-5a,b  show  the  comparison  between  experi- 
ment and  calculation.  Both  pressure  amplitudes  were  normalized 
at  the  maximum  pressure  peak.  In  both  cases,  there  are  some 
quantitative  differences,  but  good  qualitative  agreement; 
most  of  the  local  maxima  and  minima  occur  at  the  proper 
positions  and  with  the  approximately  correct  amplitudes. 

Figures  IV-6a,b  are  enlarged  graphs  of  Fig.  IV-5a,b 
for  the  region  0 <_  x <_  2xq.  From  these  graphs,  it  is 
obvious  that  the  tails  of  major  peaks  of  the  experiment 
and  theory  agreed  very  well. 
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2.  Phase 


Theoretical  calculations  were  also  done  with  the 
computer  program  "WEDGEO"  for  a distant  source.  The 
results  were  compared  with  experimental  data  (Fig.  IV- 7 , 

CASE  2 only)  with  both  phases  normalized  at  x = XQ.  In 
the  region  XQ  <_  x £ 2Xq,  these  phases  were  almost  identical, 
but  some  differences  were  observed  for  x < Xq;  the  experi- 
mentally observed  phase  changed  somewhat  more  rapidly 
than  that  of  the  theoretical  prediction,  the  difference 
increasing  as  the  apex  was  approached.  There  are  some 
possible  explanations  for  this.  The  depth  of  the  wedge 
at  x = XQ  was  4.7  mm,  and  most  differences  were  observed 
where  the  depth  was  3.0  mm  or  less  (0  x 6 cm)  . In 
this  region,  the  pressure  amplitude  was  very  small 
(Fig.  IV-2a) , so  that  measurement  of  the  phase  was  quite 
difficult.  Also,  very  close  to  the  apex  the  surface  of 
the  silicon  oil  might  be  affected  by  surface  tension  and 
not  remain  ideally  flat. 

The  incident  angles  0 in  the  upper  layer  were  approxi- 
mated as  follows;  from  Fig.  II-l,  the  angle  of  incidence 
is, 

k , 

9 = arc  cos  ) 


where 


k 


2 TT  f 
C1 
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x 2 for  which  kx,X'  have  the 


r 


Two  adjacent  points  , 
same  values  in  the  interval  -it  to  +tt  yield  a change  in 
phase  of  2 n , 


-k  ,X'  + k ,X' 
x2  2 1 


= 2 7T 


2tt 


X1  - 


2u 

Ax 


and 


0 


arc  cos 


( IV- 2) 


Figures  IV~8a,b  show  the  results  for  the  experiment 
compared  to  the  results  obtained  from  the  method  of 
images  and  from  the  simple  model.  There  is  no  turnaround 
point  (for  which  we  would  have  9 = 90°)  . In  both  cases 
the  incident  angle  equals  the  critical  angle  a little 
further  away  from  the  apex  than  x/XQ  = 1,  which  is  consis- 
tent with  the  measured  behavior  of  the  amplitude.  For 
x > XQ,  most  of  these  incident  angles  were  less  than  the 
critical  angle  0 , (This  does  not  imply  that  there  is  no 
energy  transmitted  into  the  bottom:  these  incident  angles 
result  from  the  phase-coherent  interferences  of  a large 
number  of  images,  the  fields  of  which  have  angles  of 
incidence  on  the  bottom  encompassing  a range  of  values, 
some  of  which  exceed  0 . This  shows  a significant  differ- 
ence between  the  method  of  images  and  the  simple  model  for 
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which  energy  can  not  be  transmitted  into  the  bottom  when 
the  incident  angle  is  less  than  the  critical  angle.) 

The  phase  as  predicted  by  the  method  of  images  and 
the  experimental  phase  is  nearly  linear  in  x.  The  ampli- 
tude varies  more  slowly,  but  perhaps  not  slowly  enough 
for  the  method  of  stationary  phase  to  be  applicable. 

C.  ATTEMPTED  RECONCILIATION 

The  predictions  of  the  simple  model  do  not  agree  well 
with  the  experimental  measurements.  However,  investigations 
of  the  pressure  field  in  the  bottom  would  be  facilitated 
if  a simple  analytical  expression  for  the  pressure  distri- 
bution along  the  wedge-bottom  interface  could  be  formulated. 

Consequently,  the  purpose  of  this  section  is  to  com- 
pare the  simple  model  with  the  method  of  images  for  those 
cases  where  the  source  is  very  distant.  The  amplitude  of 
the  pressure  distribution  predicted  by  the  simple  model 
has  the  form 

p.  k.  tan6  k-,  , k 

A = AQ  e ° ( IV- 3 ) 


where  the  predicted  values  are  k 
_ 3 


1,  k3  = 


1 . ’T 


and 


1 2 

^ Changing  0^,  0^*  c^,  c.,  and  b it  was  found  that 


the  pressure  amplitude  depended  only  on  the  ratios  — , 
tan  0_  °2 


Next  comparing  the  predictions  of  the  method  of  images 


with  the  simple  model  also  showed  that  the  value  k^  = 4 
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3*4* 


gave  the  proper  spatial  behavior,  to  reasonable  accuracy. 

This  was  accomplished  by  taking  In  of  both  sides  of 

A 

Eq.  (IV- 3)  and  plotting  then  In  x vs  vr  on  log- log  paper. 

o 

Obviously  the  slope  of  the  equation 

A p.  k1  tan  8 k , 

In  (In  (-— ) } = ln{(jt)  X- ( g— -)  -k^  + k4  ln(J-,) 

(IV-4) 

gives  k.. 

P1 

By  the  same  technique  the  exponents  k^  and  k2  of  — 
tan  0 p2 

and  and  the  factor  k,  were  examined  but  in  these 

3 3 

cases  the  results  for  k^,  k2,  k^  were  not  constants;  but 

p 1 C1  x ' 

functions  apparently  dependent  on  — , — , 3,  Tj-r. 

P2  C2  o 

We  can  therefore  conclude  that  the  simple  model  is 
not  generally  valid  for  arbitrary  values  of  p^,  P2 t c^, 

X * • 

c 2>  B and  ^-p.  However,  it  is  useful  to  note  that  by  assuming 

1 o 

k^  = j the  simple  model  predictions  are  brought  into  good 
agreement  with  the  experimental  measurements  for  the  values 

= 1.015  gr/cm3 

P2  = 0.995  gr/cm3 

c^  = 1173.5  m/sec 

c2  = 1503.8  m/sec 
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38.706 


3 = 0.0476  rad 


After  the  discovery  of  this  agreement  an  effort  was  made 
to  specify  the  ranges  of  variables  in  which  the  amplitude 
prediction  of  the  simple  model  is  valid  with 

3 4 

The  simple  model  requires 


x ' 1 2 

7T  <<  n COS  9 

Xo  2 c 


( IV- 5 ) 


Assuming  that  the  amplitude  distribution  is  unimportant 

X 1 A 1 

for  those  values  ^-r  such  that  j—  <<  — , this  inequality  must 

‘ o o “ 

be  achieved  within  the  restriction  (IV-5) . 


p.  tan  9 
1 c 


1 2q 
cos  9 


( IV- 6 ) 


3 << 


sin  9 ( 1 - sin  9 ) — 


( IV- 7 ) 


This  inequality  was  tested  on  the  interval  0.002  < S <_  0.005 

P 

f°r  = 60°  and  — 0.5  or  1.0  there  was  uniformly  good 

agreement  between  the  simple  model  and  the  method  of 

13 1 Pi 

images  for  — = 0.5,  but  for  — = 1.0  there  was  strong 
P2  P2 

disagreement  between  the  two  models . 


» 


p 


Therefore  the  range  of  validity  of  the  simple  model 
with  ^3  53  j appears  badly  limited.  Further,  since  the 
phase  of  the  simple  model  does  not  agree  with  the  experi- 
mental measurements  and  the  predictions  of  the  method  of 
images  (as  already  has  been  mentioned  in  part  IV-A2) , 
we  can  conclude  finally  that  this  model  is  very  weak  and 
its  application  very  restricted. 
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Kig./V-la  Experimental  pressure  distribution  CASE 


N0KMAL1ZED  DISTANCE  a/X 

Klg./v-2a  Comparison  of  pressure  amplitude  i 
(Experiment  VS.  simple  model) 


Fig./V-2b  Comparison  of  pressure  amplitude  CASK  2 
(Experiment  VS.  simple  model) 


P’ig./V-4  Theoretical  phase  and  Incident  angle  (Simple  model) 


NOKMALIZh'D  DISTANCE  -V/Xo 

Comparison  of  pressure  amplitude  GAi 
(Experiment  VS. method  of  images) 


NORMALIZED  DISTANCE  A'/Xo 

Kig./V-8a  Comparison  of  incident  angle  CASE  2 
(Experiment  VS.  method  of  images) 


: 


i 


V.  CONCLUSIONS 

Two  main  conclusions  should  be  emphasized. 

(1)  After  a careful  examination  of  the  simple  model 
it  is  found  to  fail  to  describe  adequately  the  pressure 
along  the  bottom  and  it  is  therefore  not  useful  for  future 
efforts  in  describing  the  pressure  radiated  into  the 
bottom. 

(2)  On  the  other  hand  the  method  of  images  gives  a 
good  prediction  of  pressure  amplitude  and  phase  distribu- 
tions along  the  interface  and  therefore  offers  us  hope  in 
attempting  to  predict  the  properties  of  the  beam  of  sound 
in  the  bottom. 
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APPENDIX  A 


A computer  program  for  the  calculation  of  the  pressure 
amplitude  and  phase  distribution  along  the  bottom  by  the 
method  of  images. 
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FLOW  CHART  OF  THE  SUBROUTINE  "PRES" 
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„ BBS!  ^tTYmCTlCAB13 

sssss-*-  " 


WEDGE  0 

***  MAIN  PROGRAM  *** 


**********♦*****#*****#**************##***♦**#***#***##*******.******* 
* THIS  PROGRAM  EMPLOYS  THE  METHOD  OF  IMAGES  TO  OBTAIN  THE  PRESSURE  * 


♦ AND  PHASE  DISTRIBUTION  ALONG  1 HE  BOTTOM  OF  A WEDGE-SHAPED  FLUID  * 

♦LAYERfOVERLYING  A FAST  BOTTOM.  * 

* * 

* ASSUMPTION  * 

* PLANE  WAVE  REFLECTION  COEFICIENTS.  * 


*♦#*******♦♦#♦**#*#***#*♦******#**,####*#*. *******♦**♦**************:*  Y..Y 

* NOTATION 

I/O  OUT-OUTPUT 
IN  =INPUT 
CHARACTER 

U VARIABLE 
P =F-AR  AMETER 
C -CONSTANT 


R -REAL  NUMBER 
I -INTEGER 
Z -COMPLEX  NUMBER 


SYMBOL 

MEANING 

I/O 

CHARACTER 

PA 

PRESSURE  AMPLITUDE  ON  THE  BOTTOM 

OUT 

V.R 

PH 

PRESSURE  PHASE  ON  THE  BOTTOM 

OUT 

VfR 

Cl 

SOUND  SPEED  IN  MEDIUM  1 ( UEDGE ) 

IN 

CfR 

C2 

SOUND  SPEED  IN  MEDIUM  2 (BOTTOM) 

IN 

CfR 

ROl 

DENSITY  OF  MEDIUM  1 

IN 

CfR 

R02 

DENSITY  OF  MEDIUM  2 

IN 

C f R 

BETA 

SLOPE  ANGLE  OF  BOTTOM 

IN 

CfR 

X 

SOURCE  DISTANCE  FROM  1 HE  APEX  IN  METERS 

IN 

CfR 

XX 

DISTANCE  FROM  THE  APEX  IN  METERS 

OUT 

VfR 

AT  UHICH  PA  AND  PH  ARE  CALCULATED 

F 

DRIVING  FREQUENCY 

IN 

C f l> 

H 

SOURCE  DEPTH  IN  METERS 

IN 

CfR 

ANGLE 

CRITICAL  ANGLE 

CfR 

HO 

LOWEST  POSSIBLE  MODE  CUT  OFF  DEPTH  IN  METERS 

CfR 

XO 

LOWEST  POSSIBLE  MODE  CUT  OFF  DISTANCE  FROM  THE 

CfR 

APEX  IN  METERS 

XN1 

DISTANCE  FROM  APEX  AT  UHICH  CALCULATION  STARTS 

IN 

CrR 

NORMALIZED  DY  DIVIDING  BY  XO 

XN2 

DISTANCE  FROM  AFEX  AT  WHICH  CALCULATION  STOPS 

IN 

CfR 

NORMALIZED  BY  DIVIDING  BY  XO 

distxi 

DISTANCE  XN1  IN  METERS 

C » R 

DISTX2 

DISTANCE  XN2  IN  METERS 

CfR 

DX 

RANGE  OF  CALCULATION  IN  METERS 

CfR 

N 

NUMBER  OF  POINTS  FOR  WHICH  PA  AND  PH 

IN 

C.  I 

CALCULATED 

PA  I 

3.1415 

CtR 

UN 

WAVE  NUMBER 

CfR 

C21 

SOUND  SPEED  RATIO  -C2/C1 

CfR 

RC21 

ACOUSTIC  IMPEDANCE  RATIO 

C » R 

ANGLO 

ANGLE  FORMED  BY  SOURCE t APEX  AND  BOTTOM 

CfR 

D TGTX 

TEMPORARY  FOR  XX(2)/X 

P f R 

PAMP 

TEMPORARY  FOR  PA 

PfR 

F'HAS 

TEMPORARY  FOR  PH 

PfR 
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non  n n n n n n n n n n n non  non  non  o n n n n n 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
F»OM  OOPY  PURaiSHjSD  TO  DJDC . 


************* *****  **********  ****************  **************************** 

< (PROGRAM)  ) 

DEFINE  CHARACTERS  OP  VARIABLES . PARAMETERS  AND  CONSTANTS 
IMPLICI T REAL*£l(A-H.O-Z> . IN TECER ( I -N ) 

DIMENSION  PA< 101 ) fPH( 101 ) .XX( 101 ) 

COMMON  cm  .RC21  .BETA.  UN.  PA  I .ANGLO.  XfHfF 

READ  IN  DATA 

READ < 5 » 1 00 ) C1.C2fRQIiR02fBETA 
RE AO ( 5 . SOO ) XN 1 . XN2 . X f H . F 
READ ( 5 . 600 ) N 

CALCULATE  AND  WR1 TE  CONSTANTS 

WRITE  ( 6 . 700 ) C1fC2fR01fR02fXN1fXN2fX.HfBETAfF 
ANGLC«DAF<COS  < C l /C2  > 

HO=C1/(4.000*F*DSIN(ANGLC)  ) 

XO=HO/DSIN( BETA) 

WRI TE< 6.400)  HOfXO 

TRANSFER  NORMALIZED  DISTANCE  TO  ACTUAL  DISTANCE 

DISTXl=*XNltXO 

0:STX2=XN2*XO 


CALCULATE  RANGE  DX 
DX -DISTX2-D1 S TX l 

CALCULATE  CONS  1 ANTS  FOR  SUBROUTINE 

PAI  = 4 . ODO*DA  T AN  < 1 .ODO) 

UN-2.0D0*PAI PF/C1 
C21 -C2/CI 

RC21 =C2*RQ2/ ( Cl *R01 ) 

ANOLO=BE 1 A-DARSIN ( H/X ) 

WRITE  ( 6 f 200 ) 

CALCULATE  PRESSURES  WITH  RESPECT  TO  DISTANCE  XX. 
WHERE  XX  IS  CHANGED  BY  I NCR 1MENT 

L -N+  1 

DO  101=1 fL 

XX< I > =D l STX 1 1 DFLOAT  < I- 1 ) 4DX/DFL0A  T ( N ) 

DISTX=XX( I ) /X 

CALL  P RES  (OISTXfPAMPfPHAS) 

PA ( I )»PAMP 
PH  < I )=PHAS 

WRITE  (6.J00)  XX< I ) fPA ( I ) fPH< I ) 

10  CONTINUE 

WRITE  (6.800) 

PLOT  THE  PRESSURE  AMPLITUDE  VS.  DISTANCE  XX 
AND  PRESSURE  PHASE  VS.  DISTANCE  XX 

WRITE  (6.700)  Cl .C2.R01 .K02.XN1 f XN2 f X f H f BETA f F 
CALL  DPLTP  (XXfPAfNfO) 

WRITE  (6.800) 

WRITE  (6.700)  Cl .C2.R01 .R02.XN1 .XN2.X.H. BETA.F 
CALL  DPLTP  (XX.PHfNfO) 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABU 
FROM  OOPY  FUKAISHEL'  TO  DDC 


I/O  t-ORMATS 


100  FORMAT 
200  FORMAT 
J 00  FORMAT 
•100  FORMAT 
C 

500  FORMAT 
600  FORMAT 
700  FORMAT 
C 

c 

0 

c 

c 

000  FORMAT 
S TOR 
END 


(5F15.0) 

<//;'.'X  , • HI  STANCE  ' . 1/X.  ' PRESSURE  AMP.  ' . 12X,  'PHASE  ANGLE'  //  ) 
( U0X.D15.7, 1OX.015. /» 10X.U15. 7) 

<//UOX.  LOUEST  POliS  I DLC  MORE  CUT  OFF  DEPTH' .‘7X.D15.7 
/'UOX.  'DISTANCE  FROM  APEX  ' » UOX  ,015. '///  > 

(5F15.0) 

( 14) 


( //USX . ' INPUT  DATA' . 15X, ' UFDGEO ' // 

USX. 'Cl  =' ,015.7, 5X. 'C2  ='.015.7/ 

USX.'ROl  = ' .IU5.7.5X, 'R02  ='.015.7/ 


USX.'XNl  =• ,015.7, 5X, ' XNU 
USX.'X  =•' .015.7, 5X, 'H 

UbX » ' DETA= ' » D IS . 7 , 5X , ' F 


=' ,015.7/ 
=' .015.7/ 
=' ,015.7/) 


( 1M1  ) 


C 
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***************  ...,.  aC.f  [MF-ur  CONSTANTS 


AND  F 


_ cjO'M  ,fP  TA.UN.PAI  »X»H 

UHE  L (’em  MAIN  RUUlINE.  IN  ROUTINE.  AND. 

•SKs  s'rs- — 


NOTATION 


symbol 

"pres 

p 


ANOLN 

BD 

R 

N 

II 

MM 

THETAO 

DR 


MEAN I NO 


’I/O  CHARACTER 


THE  S AM 

CHECK 
PSA  1 
REEL 
Rfc.PL  NO 
OIRET 
Z 


^RE^rmrCAUSEr BT  N TH^AIR  fM^OEtSlMAGES  ARE 
SUM  OF  PRES. < IK  "LL 

p >MAGE,  APEX  and  BOTTOM 

PARAMETER  ««S-SM8iESE<H',TH  IMAGE 
S0Trf  iSS^S  "t?R.b  ok  -oes 

^Ax"MUMrNUMBF.RSNOK  BOTTOM  BOUNCES^  reCEIVInG  POINT 

[MAGE  TO  the  RECEDING  POINT 
“no^uzed  IVtMUIOINGt^/0 

°\1^rERMEDIATE  REFLECTION^  g crITXCAL  ANGLE 

SmFTER^ REFLECTION  COEFFICIENT  B0UMCE 

Sssss;  s *»  ~~ 




P.Z 

P.Z 

P.R 

P»R 

P.Z 

P.I 

FfI 

PfI 

P.R 

PfR 

PfR 

PfR 

PfZ 

P.Z 

P.Z 

P.R 

P.Z 
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I 


Is  B£ST  quality  pkacticabu 

«»«  OOhY  JURHISHSD  TO  DOC 


c ************************************** *********** ************ ********* 

C 

C < ( FKUGRAM ) ) 

C 

SUBROUTINE  PRES  UUSTXfPAMPfPHAS) 

C 

C DEFINE  CHARACTERS  OF  VARIABLE fPARAME 1ERS  AND  CONSTANTS 
C 

IMPLICIT  REAL*H(A-HfO-Z) . INTECER(I-N) 

COMPLEX* 16  PSAI fREFLfRLFLNOfRfZ.pRESfP 
COMMON  C21 fRC21 f BETA. WNfPAIf ANGLO fXfHfF 
C 

£ INITIALIZE  COMPLEX  PRESSURE  P=0.0+J0.0 

P’PCMPLX  < 0 . ODO  f 0 . ODO ) 

C RESET  COUNTER  N»0 

c 

N=0 

10  CONTINUE 
C 

C CALCULATE  ANGLN  AND  DETERMINE  IF  ANGLE  IS  IN  THE  RANGE  OR  NOT 
C IF  ANGLE. GT.PAI  GO  TO  THE  NEXT  STEP 
C 

ANGLN=2 , ODO* I NT  < (N+l . >/2. > *BET  A+ ( - 1 ) **N* ANGLO 
IF  <ANG(,N.GE.PAI>  GO  TO  50 
C 

C CALCULATE  PARAMETER  THETAO  AND  DR 
C 

DD=DCOS  < ANGLN ) -D I S T X 

THE T AO  ~DA  T AN2 ( DSIN < ANGLN ) fDD) 

DR=DSQR  T ( 1 . ODO+DISTX**2-2 • ODO*DIS  TX*DCOS ( ANGLN ) ) 

c 

C INITIALIZE  THE  TOTAL  REFLECTION  CQEKFCIENT  R=1.0+J0.0 
C 

R=DCMPLX ( 1 .ODOfO.ODO) 

C 

C DEFINE  MM  BY  EACH  PATH  WHICH  IS  DETERMINED  BY  THE  N TH  IMAGE 
C 

MM=  I D I N T ( ANGLN/  ( 2 • QDO*DET  A ) ) 

L=MM+1 
DO  40  M=1fL 
C 

C CALCULATE  PARAMETERS  THE'TAM  AND  CHECK 
C 

THET AM= THET AO— 2 . 0D0*0FLGA  T < M-l ) *BETA 
CHECK= 1 . 0D0-C21 **2*DC0S  < THE TAM  > **2 
C 

C IDENTIFY  THE  INCIDENT  ANGLE  WHICH  IS  LESS  THAN  THE  CRITICAL 
C ANGLE  OR  NOT 
C 

IF ( CHECK. GT.O. ODO)  GO  TO  20 
C 

C CALCULATE  PARAMETER  PSAI  FOR  THE  REFLECTION  COEFFC IENT  OF  EACH 
C DOUNCE.  THERE  ARE  TWO  WAYS  WHICH  DEFEND  ON  THE  IDENTIFICATION 
C OF  A CHECK 
C 

PSAI^DCMPLX  < 0 . ODOf -DSOR  T< -CHECK ) /DS IN < THETAM  > ) 

GO  TO  30 
20  CONTINUE 

F SA I = DCMPLX  < DSC1RT  ( CHECK ) /DS IN < THETAM  > f 0 . ODO  > 

30  CONTINUE 
C 

C CALCULATE  REFL 
C 

REFLa<  RC21-PSAI ) / < RC2 1+F'SAI ) 
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TTU3  PiCTC  IS  BEST  QV1ALVTC  rKACTlCABLI 

n*va  ixm  njk»ishw»  ro  iu>c 


c 

c it  nit  POUNCt  is  nit  lasi  r hum  mt  sourci:  . mi  n mu  no-mil 
c 

niM.to.n  mh.no-m.ii. 

c 

C CM  Oil. A ft  INC  IUMI  MILtCTlIiN  CUtTl'C  l IN  f >5  I OK  I AIM  I'AIM 

C 

15  -RtKKFl. 

•*0  CONTINUE 

C 

C CMLCIIl  A 1 t I’AKAMI  r t KS  ANU  I'MS 
c win  m in  in  i is  mi  nihttiuvin  nr  rut  source 
o mm  t is  usi.n  iinl i whin  mi  soukci:  is  ii in o i ionai 
C 01  lit  K Wist  lUMI-l.O 
C 

.•'iu:m  lxio.ouo.  wn»\«I'N' 

IUM.I  "DSlNl  0.0101>0*UN*HSIN»  I HI  r Ah  ANSI  0>  */ 

C (O.OlOllOtUNnilSlNl  IKtTAM-ANUl.O)  * 

PMS-IUM r *1  in  xn  / I Hi  l > ** INI ( vNI 1 . > 
c:  1 1 .oooi  i > o iu)  mu  noi *ii-  oh 

C 

c I'Aicm.Ait  nit  mossum.  p which  is  a sum  oi  nits 
c 

p-p  it-ms 
c 

0 Str  III  lilt  NfXt  PA  1 15  OF  IMAlit 

c 

N-Ntl 

1)1)  ID  10 
SO  CONIlNUt 

0 

C CAt.l'l.’l  AIL  I'A  AND  I II.ANH  M IIIKN  PACN  TO  lilt  MAIN  KOI) T INF 

c 

I AMI' -OlWll<S  V I' ' 

I I. I AL  *P 
P 't'-PAtAl. 

f r •■*ln‘MI  I Xd'.t'OO.  I ,Ml>0> 
r I MAll  «PS  IQNl  Cl'AHS  iK'i  -Pp  ' 
t HAS-UA  I AN.'  I I I MAW  .I'M  AL  > 

m iiikn 


APPENDIX  B 


INCIDENT  ANGLE  FOR  A DISTANCE  SOURCE 

General  expression  of  an  incident  angle  is  given 
by  Eg.  (II  39-a) 


e 

nm 


0 _ - 2m8 
no 


(II  39-a) 


where  the  maximum  number  of  bounces  M is  limited  by 


4> 

M - INT 


(11-40) 


For  a distance  source 


0 

mo 


(11-48) 


A few  incident  angles  8 for  a distance  source  are  shown 
in  Table  B-l.  From  this  table  it  is  seen  that  the  reflec- 
tion coefficient  for  the  4th  image  is 


[(4)  - R(^40)R(e4i)R(t?42)  - r ( 2 ) r ( t*4Q ) 

and  by  mathematical  induction  Eqs.  I 1-50  can  be  deduced. 
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